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Abstract 

We present a combined experimental and theoretical study on the role of copolymer architecture 
in the self-assembly of binary PEO-PCL mixtures in water-THF, and show that altering the chain 
geometry and composition of the copolymers can control the form of the self-assembled struc- 
tures and lead to the formation of novel aggregates. First, using transmission electron microscopy 
and turbidity measurements, we study a mixture of sphere-forming and lamella-forming PEO-PCL 
copolymers, and show that increasing the molecular weight of the lamella-former at a constant ratio 
of its hydrophilic and hydrophobic components leads to the formation of highly-curved structures 
even at low sphere- former concentrations. This result is explained using a simple argument based 
on the effective volumes of the two sections of the diblock and is reproduced in a coarse-grained 
mean-field model: self-consistent field theory (SCFT). Using further SCFT calculations, we study 
the distribution of the two copolymer species within the individual aggregates and discuss how this 
affects the self- assembled structures. We also investigate a binary mixture of lamella- formers of dif- 
ferent molecular weights, and find that this system forms vesicles with a wall thickness intermediate 
to those of the vesicles formed by the two copolymers individually. This result is also reproduced 
using SCFT. Finally, a mixture of sphere-former and a copolymer with a large hydrophobic block 
is shown to form a range of structures, including novel elongated vesicles. 
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I. INTRODUCTION 



Amphiphiles such as block copolymers and lipids can self-assemble into many different 
structures when dissolved in solution The case of block copolymers has proved es- 

pecially interesting to researchers in recent years, for a variety of reasons [3]. First, the 
study of block copolymers is a promising route to a fundamental model of the self-assembly 
of amphiphiles in solution, since the theoretical understanding of the constituent polymer 
molecules is on a firm footing. Well-established methods such as self-consistent field the- 
ory (SCFT) 4-0] have provided considerable insight into the self-assembly of polymers, 
especially in melts (3, Isj], whilst using simple models of the individual polymer molecules. 
Second, vesicles formed from block copolymers show more promise as vehicles for drug de- 



livery |9| than similar structures formed from lipids, as the thickness and low solubility of 
their membranes means that they can 

For solutions of a single type of diblock copolymer, it is often relatively straightforward 
to understand why a given type of aggregate forms in a given system. The main factor 
that determines the shape of the structures is the architecture of the copolymer; that is, 
the size of its hydrophilic and hydrophobic blocks 12I (although other factors, such as the 
overall size of the copolymer, may also play a role [13J). If the hydrophilic component is 
large compared to the hydrophobic component, then curved aggregates such as spherical 
or cylindrical micelles form. Conversely, if the hydrophobic component is large, lamellar 



structures such as vesicles are observed [12j|. Recently, we demonstrated this behavior in 
a study of PCL-6-PEO block-copolymers with various volume fractions of the hydrophobic 
block (PCL) 14J. Here, large volume fractions /eo of the hydrophilic block (PEO) resulted 
in micelles (/eo > 0.3), lower /eo favored wormlike micelles (0.25 < /eo < 0.3) and still 
lower fractions (/eo < 0.25) led to the formation of vesicles. 

We can gain increased control over the self-assembly by mixing two types of amphiphile 
that individually form aggregates of different curvatures. Such mixtures are well known 
in cell biology, where different lipids can be sorted by segregation to regions of high and 



ow curvature 
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la] , and have also been studied in the context of lipid-detergent systems 



181 ] . More recently, they have been investigated in block copolymer solutions. For 
example, Jain and Bates [2j have studied mixtures of polyethylene oxide-polybutadiene 
(PEO-PB), and have found that blending ratio can be used to control self-assembly, and 
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furthermore that novel structures such as undulating cylinders form. They also found that 
different aggregates form depending on whether the two polymer species are mixed before 
or after their individual self-assembly j^J. 

In a recent study of a mixture of sphere-forming and lamella-forming polycaprolactone- 
co-polyethylene oxide in water-THF mixed solvents 19) , we have built on this work by con- 
trolling the quantities of water and THF to mix sphere- and lamella-forming copolymers not 
only before and after but also during their individual self-assembly Those copolymers mixed 
before self-assembly (pre-mixed) formed a sequence of aggregates of increasing curvature as 
the amount of sphere-former was increased, forming vesicles, then a mixtures of vesicles, 



rings and worms, and finally spherical mice 



been observed in lipid-detergent mixtures 



sel 



-consistent field theory 



13, 



les. This series of shape transitions has also 



18) . and has been studied theoretically using 



19M21) and models of chain packing 22J and membrane curvature 



23) . When mixed after self-assembly (post-mixed), the two species remained locally in the 
equilibrium states of the pure components, and a mixture of vesicles and spherical micelles 
was observed. The structures found when the two species were allowed partially to self- 
assemble before mixing (intermediate mixing) were more unusual, and included metastable 
paddle- and horseshoe-shaped aggregates. Using self-consistent field theory, we reproduced 
the transitions between morphologies observed in the pre-mixed system and also details of 
the aggregates such as the bulbous ends of the rods [19j]. We also gained insight into the 
complex metastable structures seen at intermediate mixing, by showing in SCFT calcula- 
tions that the segregation of the two types of copolymer can stabilize regions of different 
curvature within a single aggregate. 

In the current paper, we extend this study by varying the architectures of the copolymer 
species. We consider three specific cases: varying the length of the lamella-former in a 
blend of sphere- and lamella-formers, blending two lamella-formers of different lengths, and 
blending sphere-former with a polymer that has such a large hydrophobic block that it 
precipitates in solution if not mixed with more hydrophilic molecules. In all cases, we observe 
the quantitative and qualitative changes in the self-assembly as the copolymer architectures 



are changed. As in our previous work 



19) , we perform self-consistent field theory calculations 



in tandem with our experiments and discuss how the distribution of the two copolymer 
species within the self-assembled structures leads to the formation of the structures seen in 
the experiments. 
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The article is organised as follows. In the following section, we give details of our exper- 
imental and theoretical methods. The Results section is divided into three subsections, one 
for each of the mixtures introduced above. We then present our conclusions. 

II. METHODS 

A. Materials 

The PEO-PCL block copolymers were purchased from Advanced Polymer Materials Inc., 
Montreal and used as received. GPC analysis was also provided by Advanced Polymer Ma- 
terials Inc. and was referenced against PEO standards. Degrees of polymerization for the 
PCL block were calculated by 1 H NMR in CDCI3 by comparison to the PEO block (the 
degrees of polymerization for the monomethoxypoly(ethylene oxides) used in these polymer- 
izations are known). The molecular weight and molecular weight distributions are given in 
Table HI All other reagents with the exception of NMR solvents were purchased from Sigma 
Aldrich Company Ltd., Gillingham. Standard solvents were of spectrophotometric grade 
and inhibitor free. Deuterated NMR solvents were purchased from Euriso-top S.A., Saint- 
Aubin. All solvents were filtered before use through Pall Acrodisc PSF GHP 200 nm filters. 
For all experiments, distilled and de-ionized Millipore water (resistivity = 18.2Mfi.cm) was 
additionally filtered through Sartorius Ministart 200 nm filters directly before use. 

B. Preparation of Solutions 

Aqueous dispersions of block-copolymer aggregates were prepared by dissolving the poly- 
mer in THF to a concentration of 10 mg ml -1 . These solutions were then mixed in the volume 
ratio noted for the experiments. All the mixing ratios are thus ratios of the masses of the 
respective polymers as opposed to molar ratios. Due to the close molecular weights (Table 
H]) of the two copolymers PCL 5 kPEOik and PCL 5 kPE02k these ratios are not very different 
in this case, while for the other combinations the conversion is easily calculated. Millipore 
water was added either manually or by an Eppendorf EDOS 5222 Electronic Dispensing 
System. 125 aliquots of 20 /d of water were added at one-minute intervals. 
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Commercial sample code 


Sample formula 


M w * 


M w /M n 


/eo b 


Morphology c 


R h (nm) d 


PCL 10k PEO 2k 


PEO45 - 6-PCL101 


17300 


1.36 


0.15 


V 


220 


PCL 5k PE0 lk 


PE0 23 - &-PCL 47 


7100 


1.15 


0.17 


V 


170 


PCL 5k PE0 2k 


PEO45 - &-PCL43 


7800 


1.16 


0.3 


c,s 


30 


PCL5 k PE055Q 


PEOi 2 - 6-PCL 56 


6380 


1.15 


0.08 


p 


n/a 



TABLE I: Polymers used in this study, with their properties from light scattering. 

Notes: (a) total molecular weight from GPC (PEO standards) (b) volume fraction of the EO block 
calculated from the melt densities of the two blocks (c) morphology determined by light scattering 
and cryo-TEM (S=spherical micelles, C = worm-like micelle, V = vesicle, P = precipitate; in the 
cases of mixed morphologies the majority component is written first), (d) hydrodynamic radius 
from dynamic light scattering (DLS) after dialysis. 

C. Turbidity Measurements 

We performed turbidity measurements during the preparation of the samples using an 
adapted Perkin Elmer UV/Vis Lambda 40 Spectrometer. A wavelength of 600 nm was used 
with a slit width of 2 nm. Stirring was performed using a standard magnetic stirrer/hotplate 
placed under the spectrometer. The polymer was dissolved in THF (1 mL, 10 mg mL _1 ) and 
a zero reading was taken (transmittance, T = 100%). Millipore water was then added either 
in 10 fil aliquots every 30 s using an Eppendorf EDOS 5222 Electronic Dispensing System 
and a turbidity reading was taken after each addition. 



D. Cryo-TEM 

Samples for thin-film cryo-TEM were loaded onto plasma-treated (30 seconds) holey- 
carbon grids and prepared using a GATAN cryo-plunge into liquid ethane and then trans- 
ferred using a GATAN 626 cryo-transfer system. Samples were examined using a JEOL 2100 
TEM operating at 200 kV. Images were obtained using a Bioscan or a GATAN Ultrascan 
4k camera and analyzed by GATAN Digital Micrograph version 1.71.38. During our previ- 



ous investigations 



14j we observed that imaging of the self-assembled structures (especially 



vesicles) is greatly improved in samples containing ca. 30% THF compared to samples in 
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pure aqueous solution. Similar structures were observed in both solvent conditions, which 
indicates that at THF fractions of 30% and below, the mobility of the block-copolymers is 
too restricted to allow for further growth of the aggregates. However, we found that, over a 
timescale of a few months, internal rearrangements occurred that evened out local variations 
in surface curvature transforming the more complex metastable aggregates into vesicles or 
nested onion-like structures. In order to focus on the initial metastable structures that are 
formed (prior to any slow internal rearrangement processes) and to obtain as high quality 
imaging as possible, we prepared the solutions for cryo-TEM in aqueous solutions containing 
28% THF and imaged these samples within a maximum time of two weeks. 



E. Self-consistent field theory 

To further our understanding of the role of polymer architecture on self-assembly, we 
performed self-consistent field theory (SCFT) calculations 4j on a model system of two 
species of AB diblock copolymers (lamella- and sphere-forming respectively) blended with A 
homopolymer 'solvent'. SCFT is a coarse-grained mean-field theory in which the individual 
polymer molecules are modeled by random walks and composition fluctuations are neglected. 



For sufficiently long polymers [24|, these approximations prove extremely effective and the 
predictions of the theory are very accurate for a wide range of systems |6|]. We use a 
simple implementation of the theory where the interactions of the polymers are included 
by imposing incompressibility and introducing a contact potential between the A and B 
monomers 6|. The fine details of the polymer molecules are not taken into account, so that 
monomers of all species are taken to have the same length a and volume 1/po- 

From a technical point of view, a self-consistent field theory calculation consists of solving 
a series of differential (diffusion) equations to calculate the density profiles of the various 
polymer species. An initial guess for the profiles is made, which has the approximate form 
of the structure that we wish to study. The density profiles are then recalculated until a set 
of equations reflecting the physical properties of the system (such as its incompressibility) 
jf], [25I is satisfied. The SCFT differential equations are solved using a finite-difference 
method (26 1 with a spatial step size of 0.04 aiV 1 / 2 , where N is the number of monomers in 
the sphere-forming species, and reflecting boundary conditions are imposed at the origin 
and edges of the box. Full technical details of our calculation can be found in a recent 
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publication [21 ]. 

In the current paper, we use SCFT calculations in two slightly different ways. First, 
we perform effectively one- dimensional calculations on spherical micelles and infinite cylin- 
ders and bilayers to calculate simple phase diagrams as a function of sphere-former volume 
fraction and reproduce the basic phenomenology observed in the experiments. These cal- 
culations proceed as follows [27|, |28j. To begin, we calculate the free-energy density of a 
box containing a single spherical, cylindrical or planar aggregate surrounded by solvent. 
The shape of this box is set by the symmetry of the aggregate; for example, a spherical 
micelle is formed at the center of a spherical box. The calculation is therefore effectively 
one-dimensional. The volume V of this simulation box is then varied, keeping the volume 
fraction of copolymer constant, until the box size with the minimum free-energy density is 
found. Provided the system is dilute, so that each aggregate is surrounded by a large volume 
of solvent and the aggregates do not interact with each other, this provides a simple model 
of a larger system (of fixed volume and fixed copolymer volume fraction) containing many 
aggregates. The reason for this is that such a system minimizes its free energy by varying 
the number of aggregates and hence the volume ('box size') occupied by each. Although 
computationally inexpensive, this approach yields accurate information on micelle shape 
transitions and its results agree well with experiment [28] . 

We also carry out more detailed calculations on the rod and ring structures seen in the 
experiment. We have two aims here. First, we wish to show that these more complex struc- 
tures can be reproduced in detail in our calculations. Second, we will study the distribution 
of the sphere-forming and lamella-forming copolymers within the aggregates. Since both 
rods and rings have cylindrical symmetry, we will perform our (effectively two-dimensional) 
calculations in a cylindrical box. We note that it is not possible to include information on 
the distribution of the sphere- and lamella-formers within the micelles in the initial guess 
for the SCFT calculations. Any segregation of the two species will therefore arise naturally 
from the theory and need not be artificially introduced [21]. 

Although this model is relatively simple, we found in our previous work on the self- 



assembly of binary PEO-PCL mixtures 



19] that it contains enough detail to yield informa- 



tion on the structures formed in such systems and on the distribution of the two polymer 
species within these aggregates. In this earlier paper, we focused on a mixture of lamella- 
forming PCL 5 kPEOik and sphere-forming PCL 5 kPE02k- We modeled the lamella-former by 
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Polymer 


Sample code 


SCFT monomers 


SCFT N A /N B 


lamella-former (long) 


PCL 10k PEO 2k 


3N/4 


1 


lamella-former (short) 


PCL 5k PE0 lk 


N/2 


1 


sphere-former 


PCL 5k PE0 2k 


N 


3 


solvent 


n/a 


N/2 


n/a 



TABLE II: Parameters of the polymers used in SCFT calculations. For each SCFT poly- 
mer, we list the corresponding polymer in the experiments, the number of (SCFT) monomers and 
the volume ratio of the hydrophilic and hydrophobic blocks. 



a symmetric AB diblock with equal numbers of monomers Aa and Ab in its hydrophilic (A) 
and hydrophobic (B) sections. For simplicity, the homopolymer was taken to have the same 
length as the lamella-former. In line with the experiments, the sphere-former contained 
the same number of hydrophobic monomers as the lamella-former but a larger number of 
hydrophilic monomers, so that = 3 Ah- The x parameter setting the strength of the 
interaction between the A and B species was set to x = 30/A^, where A" is the total number 
of monomers in the sphere-forming species. Our aim here was not to match the experimen- 
tal polymer parameters exactly, but to reproduce the basic phenomenology of the system 
(sphere- and lamella-forming species, matched hydrophobic blocks) as simply as possible. 
We take a similar approach in the current paper. To study the effect of lamella-former 
length on the blend of sphere- and lamella-former, and to observe the result of blending two 
different lamella-formers, we introduce the larger lamella-forming copolymer PCL 10k PEO2 k . 
We model this new molecule by an SCFT polymer with 3N/4 monomers, while keeping 
Aa = Ab. The number of monomers 3N/4 is chosen since increasing the size of the sym- 
metric copolymer too much will lead to its forming micelles rather than bilayers Ij], llo] ]. 
For the sake of clarity, we summarize the SCFT polymer parameters introduced above in 
Table [III 
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III. RESULTS 



A. Blends of lamella- and sphere-formers 

To begin, we consider the structures formed when various concentrations of the sphere- 
forming copolymer PCLskPEC^k are added to the long lamella-former PCLiokPE02k- These 
molecules are chosen to isolate the effect of molecular weight on the shape transitions: they 
contain more monomers than the PCLskPEOik copolymers studied previously, but have the 
same ratio of hydrophobic to hydrophilic blocks. We first turn our attention to the turbidity 
traces of this system. The features in a turbidity trace can be directly linked to the points 
where the transitions between spherical micelles, wormlike micelles and vesicles occur jl^ . 
Specifically, in clear solutions spherical micelles or short worms dominate and no or very 



few vesicles are present. Conversely, high turbidity of the so 



ution at high water content 



indicates the presence of larger aggregates such as vesicles 14(. From these measurements 
(Figure ED it can be seen that a mixture of 10% PCL 5k PE0 2 k with 90% PCL 10k PEO 2 k does 
not behave very differently from a system of pure lamella-former, since the traces for the 
two systems are very similar. However, a sharp change in the optical transmission is seen 
on addition of 20% PCLskPEC^k- Here, although the concentration of sphere-former is still 
relatively low, the trace resembles that of pure PCL 5 kPE02k much more closely than that 
of the pure lamella-former PCL 10 kPEO2k and does not show the strong turbidity in the 
water-rich area linked with the presence of larger aggregates [l4] . We note in passing that 
the sharp dip in the optical transmission between 20 and 30% water is most probably due 
to a miscibility gap in the PEO-THF-water phase space Q, Q, [|(3] and is not associated 
with a transition in the shape of the aggregates. 

These results are in contrast to those obtained for mixtures of PCL 5 kPEOik and 

n 

PCL 5k PE0 2 k in our previous publication [19j. There, vesicles were formed up to approx- 
imately 30% sphere-former, whereas in the current system the transition from vesicles to 
micelles occurs between 10% and 20% PCLskPEC^k- Increasing the length of the lamella- 
forming copolymer is therefore seen to favor the formation of more curved structures. 

To gain more detailed insights into the system, we now consider cryo-TEM images taken 
at a range of sphere-former concentrations. In mixes with 90 % PCLiokPEC^k and 10 % 
PCL 5k PE0 2 k, this technique reveals the presence of a large variety of structures, as can be 
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FIG. 1: Turbidity traces for the self-assembly of block copolymer mixtures of PCLiokPEO^k and 
PCL5i c PE02k- The optical transmission (in %) at 600 nm is plotted against the water concentration 
in the solvent. 

seen in Figure [2j First, we note that a significant number of vesicles is still present, account- 
ing for the high turbidity of this mixture seen in Figure [TJ In addition, wormlike micelles 
and rings (end-to-end joined worms) can be seen in the left image. The wormlike micelles 
here often form branched network structures showing multiple three-way connections. The 
individual branches of the network also tend to be rather short, terminating in enlarged 

n 

end-caps. These images closely resemble some of those shown by Jain and Bates [2( for 



PEO-6-PB block-copolymer mixtures as well as those of Chen et al. [3l| for PS-6-PAA. In 
different regions of the same TEM grid, very unusual vesicles could also be seen that on 
drying deformed to create a space-filling tessellated pattern (Figure El right image). The 
majority of these structures clearly show the dark outer ring of the vesicle wall confirming 
that they are indeed closed vesicles. However, on the top right there are some structures 
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FIG. 2: Cryo-TEM of a mixture of 90% PCLio k PE0 2 k and 10% PCL 5k PE0 2 k self-assembled by 
solvent exchange from THF. The images were taken in an aqueous solution with 28% THF; the 
scalebars are 500 nm. 



that do not have this pronounced darker rim and may therefore be unwrapped bilayer sheets. 
Such structures have been proposed as intermediate stages in the self-assembly of vesicles 



32 



33|. 



At a mixing ratio of 80% PCL 10k PEO 2k and 20% PCL 5k PE0 2k (Figure 0j), the TEM 
shows a mix of spherical micelles, short worms and toroidal rings. The vesicles and sheet- 
like structures shown in Figure [2] no longer appear at this concentration. This is in line with 
the turbidity trace results: this solution is clear at high water concentrations, consistent 

n 

with the presence of small micelles such as rings and short rods |14J |. 

We now present our self-consistent field theory calculations on our simple model of our 
experimental systems. For each set of calculations, we first consider the PCL 5k PEOi k - 



PCL 5k PE0 2k mixture studied in our previous paper [19|, modeled, as described in the 
methods section, by a polymer blend including relatively short symmetric lamella-formers 
containing N/2 monomers and sphere- formers containing N monomers. Next, we move on 
to our model of the PCLio k PE0 2k -PCL5 k PE0 2k system of the current paper, in which the 
SCFT lamella-formers are still symmetric but now contain 3iV/4 monomers, and look at 
how the self-assembly is altered by the change in polymer architecture. 
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FIG. 3: Cryo-TEM of a mixture of 80% PCLio k PE0 2 k and 20% PCL 5k PE0 2 k self-assembled by 
solvent exchange from THF. The images were taken in an aqueous solution with 28% THF; the 
scalebars are 500 nm. 



To begin, we calculate the free-energy densities of ideal spheres, infinite cylinders and 
infinite bilayers using the method of variable subsystem size described above, and determine 
how these vary as the volume fraction of sphere-former is increased. To ensure that the 
system is relatively dilute and that aggregates are surrounded by a large volume of solvent, 
we fix the overall volume fraction of copolymer to 8 %. All free energy densities F/V are 
plotted with respect to that of the homogeneous solution with the same composition F h /V; 
that is, we plot the quantity / = F/V — F h /V. Since the free-energy densities /, of the 
three shapes of aggregate are quite close together, they are plotted normalized with respect 
to the magnitude of the free-energy density |/c| of the cylindrical micelle to show the shape 
transitions clearly. The cylinder free-energy density then appears as a horizontal line at 
/ = — 1, and is approached from above and below by the sphere and lamella free-energy 
densities as the sphere-former concentration increases. 

Figure H] shows the free-energy densities of spherical, cylindrical and planar aggregates 
plotted against (j)'/<p, where 0' is the volume fraction of sphere-former and <fi is the total 
volume fraction of copolymer (sphere formers plus lamella-formers). Panel a shows the 



results for the system with the shorter lamella- former of N/2 monomers In this system, 
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the lamella-former has the lowest free energy at lower sphere-former volume fractions. At 
around <f)'/4> — 35%, the lamellar and cylindrical free energies cross, and the cylinder has 
the lowest free energy until <p'/<f) ~ 40%, when the spherical micelle finally becomes most 
energetically favorable. This reproduces the series of transitions from vesicles to cylindrical 
micelles (worms and rings) to spherical micelles seen in our experiments [3], although the 
values of (f>'/4> at which the transitions occur are slightly shifted, as would be expected in 



view of the simplicity of our model. In these TEM images 19], vesicles were seen at 5% 
and 25% sphere-forming copolymer and spherical micelles are seen at 75% sphere-former, 
in line with our calculations. However, a mixture of worms, rings and vesicles is seen at 
50%, where our calculations predict spherical micelles. It is interesting to note that both 
our calculations and experiments demonstrate that a blend of sphere-forming and lamella- 
forming amphiphiles can form cylindrical micelles, even though this structure is favored by 
neither of these molecules individually [21 1 . 

In panel b of Figure HI we plot the corresponding free-energy densities for the same 
system but with the degree of polymerization of the symmetric lamella former increased 
to 3N/4 monomers: a simple representation of the longer PCLiokPE02k molecules of the 
current system. Again, we find the same sequence of morphologies as the sphere-former 
concentration is increased. However, as in the experiments, the vesicle-cylinder and cylinder- 
sphere transitions are shifted to lower sphere-former volume fractions with respect to the 
previous system l^]. Specifically, the former transition now occurs at approximately <fi' /4> = 
28%, whilst the latter is moved to around (j)'/4> = 35%. Again, due to the simplicity 
of our theoretical approach, these transitions are not perfectly aligned with those in the 
experimental system. For example, the TEM images taken at a mixing ratio of 20% (Figure 
|3]) show a mixture of rings, short rods and spherical micelles, whereas our theory predicts that 
this mixture will form vesicles, or, equivalently, that the region of coexistence of morphologies 
would be expected at a mixing ratio closer to 30%. 

We believe that the shift in the shape transitions to lower sphere- former volume frac- 
tions for increasing chain length of the lamella-former is primarily due to the greater shape 
asymmetry of the lamella former when its length is increased whilst keeping the ratio of the 
hydrophilic and hydrophobic blocks constant. This can be seen from the following heuris- 
tic model of the shape asymmetry of diblock copolymers. We assume that the hydrophilic 
block A of the lamella- former is swollen by solvent and so has an end-to-end distance 34 1 
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FIG. 4: Free-energy densities as a function of mixing ratio for (a) a short symmetric lamella- 
former with N/2 monomers blended with a iV-monomer sphere-former and (b) a longer symmetric 
lamella- former with 3iV/4 monomers blended with a iV-monomer sphere- former. The free-energy 
densities for spherical micelles, cylindrical micelles and flat bilayers are plotted with stars, circles 
and squares respectively. The transitions between the different micelle morphologies are indicated 
by vertical dashed lines. 

of i?A ~ iV A ^ 5 and effective volume Va ~ R\ ~ Nj/ 5 - In contrast, the hydrophobic B-blocks 
are taken to be in a collapsed, brush-like state 35j and so have effective volume Vb ~ N-q. 
Since we keep the ratio of the hydrophilic and hydrophobic blocks constant, we can also 
write Va ~ iV tot and Vb ~ N tot , where iVtot = + N B is the total length of the lamella- 
former. Denning the shape asymmetry of the lamella-former to be e = Va/Vb, we have 



e ~ 



Ni 



4/5 



tot 



Therefore, as we increase N tot at constant hydrophilic to hydrophobic ratio, e will 



also increase and the lamella-former will become more asymmetric. In consequence, we ex- 
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pect it to have an greater tendency to form curved structures, with the swollen A-blocks on 
the outside of the curved surface. This shift towards more curved aggregates as the overall 
copolymer length is increased has indeed been seen in experiments on symmetric diblock 
copolymers in solution 13]. 

In the above calculations, we have studied only infinite cylinders. However, the cylindrical 
micelles seen in the experiments have the form of rings or short, bulbous-ended rods. We 



now use effectively 2d SCFT calculations in cylindrical polar coordinates [19|, |2l| to study 
these structures in more detail. First, in Figure |5j we demonstrate that these structures can 
indeed be reproduced within our model's parameter space. We focus on a system of lamella 
former with 3N/4 monomers mixed with sphere- former in a ratio of 2 to 1: a blending ratio 
where cylindrical micelles have the lowest free energy in our calculations (Figure H]). To show 
the form of the structures, the sum of the hydrophilic block densities of the two species is 
plotted. Panel a of Figure [5] shows the bulbous end of a rod 3fJ, preceded by a thinner neck 
region of negative curvature [37|, [38(. We find also that ring-like structures exist as local 
solutions to SCFT, and plot a cross-section through one of these aggregates in panel b. 
We now turn our attention to the distribution of the two copolymer species within the 
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21] r)(r), which we 



rod structures. To this end, we introduce an enhancement factor 
define as _ 

,(r) = ££2|L (!) 

M r ) 0B2 

Here, 0b ( r ) is the local volume fraction of lamella-former hydrophobic blocks and 0B2( r ) 
is the corresponding quantity for the sphere-former hydrophobic blocks. The mean volume 
fractions of the two species are denoted by B and B2 . The enhancement factor tells us how 
much the concentration of sphere-former is enhanced with respect to that of the lamella- 
former at a given point in the system. We define r/(r) such that it is normalized with 
respect to the mean volume fractions of the two core species, so that values greater than one 
represent enhancement of the sphere-former concentration and values less than one represent 
depletion. The enhancement factor is plotted only within the core of the micelle, defined 
as the region where the total density of the hydrophobic species 0b ( r ) + 0B2( r ) is greater 
than that of the hydrophilic species 0a ( r ) + 0A2( r )- To locate this region accurately, we 
apply simple bilinear interpolation to our SCFT data to make the grid finer before plotting 
i](r). In Figure M (a), we show the enhancement factor for the blend of sphere- former and 
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FIG. 5: Bulbous-ended rod (a) and ring (b) structures in a system with a long symmetric lamella- 
former with 3N/4 monomers blended with a N- monomer sphere- former in a ratio of 2 to 1. The 
total density of the hydrophilic A-blocks is plotted in cylindrical polar coordinates 



shorter lamella-former (N/2 monomers) 19j. In this system, the sphere- formers segregate 
to the end of the rod 19], which is the most strongly curved part of the aggregate and, in 
fact, closely resembles a spherical micelle. Since the hydrophobic blocks of the two species 
are matched and so can both reach the center of the micelle, the ratio of the concentrations 
of the two species varies rather little in the main cylindrical body of the aggregate away 
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from the endcaps. This is in contrast to the system shown in panel b of Figure |6j where 
a longer lamella- former of 3N/4 monomers is used. Here, the sphere- formers have shorter 
hydrophobic blocks than the lamella-formers, and no longer reach the central axis of the 
cylindrical micelle. This means that the sphere-former concentration is enhanced on the 
surface of the body of the cylinder as well as in the endcap, and is strongly depleted in the 
center of the rod. We note also that this mismatch between the hydrophobic blocks means 
that segregation is a stronger effect than in the previous system, with the range of values 
taken by r/(r) significantly increased. 

To study the distribution of the sphere- and lamella-formers within the micelles in more 
detail, we now plot cuts through the density profiles of the various species both along and 
perpendicular to the axis of the cylinder for the two systems shown in Figure The plot 
for the system with shorter lamella-formers (Figure [7J confirms two points made above. 
First, it can be seen from the cut along the rod axis shown in the main panel of Figure 
[7J that the segregation of the sphere-formers to the endcap is relatively weak and that the 
peak in sphere-former concentration here is quite small. Furthermore, from the inset to 
Figure [7J which shows a cut perpendicular to the rod axis at the center of the micelle, 
we see that there is no segregation of the two polymer species in the main body of the 
aggregate. In corresponding plots for the system with longer lamella-formers (Figure IS]), we 
see strong segregation of the sphere-formers to both the endcaps and surface of the cylinder. 
In addition, plotting the data in this way also shows us that, despite the new effect of the 
enhancement of sphere-former concentration on the surface of the cylinder, the strongest 
segregation is still to the most highly-curved region: the endcaps. To see this, note that the 
peak in sphere-former concentration at the end of the rod (main panel of Figure |8]) is higher 
and more pronounced than that at the surface of the cylindrical section of the rod (inset to 
Figure [8]). We note also that the sphere-formers have a dip in concentration just before the 
peak at the cylinder endcap. This feature corresponds to the negatively-curved neck of the 



micelle 38[ visible in Figure EJ Since the sphere- formers naturally prefer positive curvature, 
they migrate away from this region. 

We believe that the stronger segregation of the sphere-formers in the mixture with the 
longer lamella-former is due to the entropic elasticity of the core blocks. Specifically, if the 
sphere-formers were homogeneously mixed with the longer lamella-formers, the hydrophobic 
cores would have to be strongly stretched, restricting the number of configurations they can 
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FIG. 6: Enhancement factor r/(r) plotted within the cores of micelles formed from (a) a 2:1 blend of 
3iV/4-monomer lamella-former and iV-monomer sphere-former and (b) a 2:1 blend of iV/2-monomer 
lamella-former and iV-monomer sphere- former. Dark areas show regions where the concentration 
of the sphere- former is enhanced with respect to the lamella- former, lighter areas show regions 
where it is depleted. 

access and so leading to a loss of entropy. In consequence, the sphere-formers move to 
the endcaps, leading to the formation of short cylinders or network structures with many 
bulbous endcaps such as those seen in Figure El 
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FIG. 7: Cuts through the density profiles of the rod structure in a 2:1 blend of iV/2-monomer 
lamella- former and iV- monomer sphere- former. The hydrophobic and hydrophilic block density 
profiles are shown with full and dashed lines respectively, and the sphere-formers are plotted with 
thicker lines. The solvent is plotted with a dotted line. The main panel shows a cut along the main 
axis of the rod and the inset shows a cut perpendicular to this axis at the center of the rod. 

The larger number of Y-junctions 3^ in the PCLiokPEC^k system compared to the 
previous PCL 5 kPEOik system [19] may be explained in a similar way jlj: in the system 
containing the shorter lamella-former PCL 5 kPEOik, the relatively short core blocks can 
adopt a smaller number of conformations and so are less able to pack into a more complex 
structure such as a Y-iunction. The increased number of branched structures at higher 

n n n 

molecular weights is in agreement with the work of several other groups [II, 13JJ, |40|, 
single- component systems. 



4l| on 
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FIG. 8: Cuts through the density profiles of the rod structure in a 2:1 blend of 3iV/4-monomer 
lamella- former and iV-monomer sphere- former. The hydrophobic and hydrophilic block density 
profiles are shown with full and dashed lines respectively, and the sphere-formers are plotted with 
thicker lines. The solvent is plotted with a dotted line. The main panel shows a cut along the main 
axis of the rod and the inset shows a cut perpendicular to this axis at the center of the rod. 

B. Blends of two lamella-formers 



Having studied the effect of blending two copolymers that individually form different 
structures, we now investigate a system where the two species are lamella-formers but of 
different lengths. The two copolymers considered are the shorter PCLskPEOik [13] and the 



longer PCL 10 kPEO2k molecule introduced above. As in our previous publication 19], we 
mixed the two polymers both before and after their individual self-assembly. 

Figure [9] shows representative cryo-TEM images of the vesicles formed by PCL 5 kPEOik 
(panel a) and PCL 10 kPEO2k (panel b) on their own. In the second row of the figure the 
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FIG. 9: Cryo-TEM images of the vesicles in 28 % THF in water. Panels a and b show the vesicles 
obtained from the self-assembly of PCL 5 kPEOik and PCL 10 kPEO2k, respectively. Panel c shows 
the structures obtained when solutions of PCLskPEOuc and PCLiokPEO^k i n pure THF are mixed 
in equal parts prior to self-assembly. Panel d shows the same composition when the two solutions 
are mixed after self-assembly at a THF content of 28%. The image was taken 1 week after mixing. 
Two distinct populations of vesicles with different wall thickness can be seen. All images are at 
the same magnification. The scalebars in a and d are 200 nm and 100 nm in panels b and c. 
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FIG. 10: Density profiles of an infinite lamellar structure in (a) an iV/2-monomer lamella- former, 
(b) a 1:1 blend of 3 N/4- monomer lamella-former and iV/2-monomer lamella-former and (c) a 3N/4- 
monomer lamella- former. The hydrophobic and hydrophilic block density profiles are shown with 
full and dashed lines respectively, and the 3N/4 lamella- formers are plotted with thicker lines. The 
solvent is plotted with a dotted line, the total density of copolymers is plotted with a dot-dashed 
line, and the point at which these two densities are equal is marked with a thick vertical dashed 
line. 

vesicles are shown that result when the two block-copolymers are mixed before assembly 
(in pure THF solution) (panel c) and after assembly i.e. after dilution to 28% THF (panel 
d). This latter case corresponds to mixing the solutions shown in panels a and b of Figure 
[9] As the chain lengths of the two copolymers differ by a factor of two, the resulting 
vesicles have different wall thickness (16 nm and 25 nm respectively). In the image d of the 
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sample mixed post-assembly, the wall thicknesses of the individual vesicles remain unchanged 
and two different populations can be clearly distinguished, while in the sample mixed pre- 



assemb 
before 



van intermediate wall thickness of ca. 20 nm is found. This result indicates that, as 



19], the exchange of material between the different structures is suppressed at THF 



concentrations below 28%. 

We now wish to see whether the result of the blending of two types of lamella-former 
leading to the formation of vesicles of intermediate wall thickness can be reproduced in our 
simple SCFT model. In panel a of Figure [T0| we plot the density profiles of the hydrophobic 
B-blocks, hydrophilic A blocks and solvent of a bilayer formed of relatively short lamella- 
formers of N/2 monomers calculated, as above, from the method of variable subsystem size. 
We also show the total density profile of the copolymers (A blocks plus B blocks) and the 
bilayer thickness, defined as the distance from the origin at which the densities of copolymer 
and solvent are equal) . Panel c shows the corresponding profiles for the longer lamella former 
with 3N/4 monomers. As expected from the respective lengths of the molecules, this latter 
polymer forms thicker bilayers than those shown in panel a. In the central panel b, we show 
the density profiles for a bilayer formed of an equal-parts mixture of the two species. In 
excellent agreement with the experiments, this mixed bilayer has a thickness approximately 
halfway between those of the pure bilayers. By comparing the density profiles in panels a 
and b we see that the addition of the longer species (in panel b) causes the core blocks of 
the shorter species to stretch outwards from their equilibrium state in a single-component 
structure (panel a). 



C. Strongly hydrophobic copolymer mixed with micelle former 

Finally, we investigate an extreme case of two strongly mismatched copolymers. The first 
of these, PCLskPEOsso, is so hydrophobic that, in isolation, it fails to assemble into vesicles 
and precipitates instead. This polymer was mixed with the micelle-former considered above, 
PCL 5 kPE02k- The two species were mixed before their individual self-assembly in a mass 
ratio of 60% PCL 5k PE0 55 o to 40% PCL 5k PE0 2k - When self-assembly was triggered by the 
addition of water, the solution became turbid, indicating the presence of large aggregates. 
This was confirmed by cryo-TEM (Figure fTTj) . which revealed the presence of a large fraction 
of vesicles with smaller populations of wormlike and spherical micelles. Remarkably, many 
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FIG. 11: Cryo-TEM of the self-assembled structures obtained for a mixture of PCL 5 kPE055o and 
PCL5kPE02k at a mass-based mixing ratio of 3 to 2 parts. The image was taken from a solution 
containing 28% THF and the scalebar is 500 nm. 

of these vesicles are not spherical but have a novel elongated shape, which could be due to 
segregation of the two copolymer species within the individual aggregates. This is likely to 
be an especially strong effect in the current system, as the solubilities of the two copolymer 
species are very different and the PCL 5 kPEOs5o chains start to aggregate at much lower water 
fractions than the more hydrophilic PCL 5k PE0 2 k- This produces a phenomenon analogous 
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to the intermediate mixing discussed in our earlier work 19[ : the strongly hydrophobic chains 
are partially self-assembled at the time they encounter the sphere-formers. In consequence, 
the two species mix less efficiently and regions of different curvature can coexist within the 
same aggregate [3]. The negative curvature regions around the centers of the elongated 
vesicles are likely to contain higher concentrations of the strongly hydrophobic copolymer, 
whereas the curved ends of these structures will probably contain more of the micelle-former. 
Due to the small difference in curvature between the different regions of the aggregate, testing 
this hypothesis is unfortunately beyond the scope of our current SCFT methods. 



IV. CONCLUSION 



In this paper, we have used a combination of experiment and theory to show that varying 
the chain geometry of binary copolymer mixtures in solution can give us precise control 
over the form of the self-assembled aggregates and lead to the formation of new structures. 
We investigated three distinct situations. First, we studied a mixture of sphere-forming 
and lamella-forming copolymers, and found in both electron microscopy experiments and 
coarse-grained mean-field theory that increasing the chain length of the lamella-former whilst 
keeping the ratio of its hydrophilic and hydrophobic components constant leads to the for- 
mation of highly-curved structures at lower sphere-former volume fractions. We presented 
an explanation of this behavior in terms of the volume asymmetry of the two sections of the 
diblock. Using more detailed SCFT calculations, we found the rings and bulbous-ended rods 
seen in the experiments, and observed a strong segregation of the sphere-forming copoly- 
mers to the curved ends of the cylindrical micelles. We explained this effect by suggesting 
that sphere-forming copolymers would pay a large free energy penalty if they were dispersed 
evenly through the aggregates, as their relatively short core blocks would need to be strongly 
stretched to fit in with those of the lamella-former. In consequence, the sphere-formers tend 
to de-mix from the lamella-formers, leading to the formation of cylinders with highly-curved 
endcaps. This segregation between species may be accentuated by other effects such as 
enthalpy of crystallization. We also studied a mixture of two lamella-forming copolymers 
of different molecular weights. In both experiment and theory, we found that the bilayers 
formed in this system have a wall thickness that is in between those observed in systems 
containing only one of the two types of lamella-former. Finally, a mixture of sphere-forming 
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copolymer and a strongly hydrophobic copolymer that precipitates in isolation was shown 
to form a range of structures, including novel elongated vesicles. 

The results presented above demonstrate the power of a combined experimental and the- 
oretical approach to the investigation and design of self-assembling block copolymers. Self- 
consistent field theory can map out the broad phase diagram of block copolymer mixtures 
and suggest experimental parameter spaces to search for new morphologies. Furthermore, 
it yields insights into the aggregates observed in the experiments, reproducing details of the 
structures and the distribution of the different polymer species within these. The wealth of 
morphologies observed in our work highlights the fine balance of forces governing the self- 
assembly behavior of block-copolymer systems. Further investigation of the different factors 
could open up a new zoo of self-assembled aggregates with the distribution and magnitude of 
local curvature differences as additional design parameters, and several avenues for extension 
of this work suggest themselves. In particular, the precise role of the sphere-former archi- 
tecture could be investigated, with the aim of producing structures of a specified curvature. 
The various components could also be mixed at different stages in their self-assembly [3], 
to access further new aggregates and to gain insight into the intermediate steps in micelle 
and vesicle formation. On the theoretical side, our SCFT calculations could be extended to 
investigate the favorability of the more complex structures, particularly the Y-junctions, in 
different copolymer mixtures. 
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